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An airborne differential absorption lidar (DIAL) system has been developed at the NASA Langley 
Research Center for remote measurements of atmospheric water vapor f H 2 0 I and aerosols. A solid-state 
alexandrite laser with a 1-pm linewidth and >99.85% spectral purity was used as the on-line 
transmitter. Solid-state avalanche photodiode detector technology has replaced photomultiplier tubes in 
the receiver system, providing an average increase by a factor of 1.5-2. 5 in the signal-to-noise ratio of the 
H 2 O measurement. By incorporating advanced diagnostic and data-acquisition instrumentation into 
other subsystems, we achieved additional improvements in system operational reliability and measure- 
ment accuracy. Laboratory spectroscopic measurements of H 2 0 absorption-line parameters were 
performed to reduce the uncertainties in our knowledge of the absorption cross sections. Line-center 
H^O absorption cross sections were determined, with errors of 3-6%, for more than 120 lines in the 
720-nm region. Flight tests of the system were conducted during 1989-1991 on the NASA Wallops 
Flight Facility Electra aircraft, and extensive intercomparison measurements were performed with 
dew-point hygrometers and H 2 O radiosondes. The H 2 0 distributions measured with the DIAL system 
differed by < 10% from the profiles determined with the in situ probes in a variety of atmospheric conditions. 
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Introduction 

Water vapor (H 2 0) is the most variable of the major 
molecular constituents of the atmosphere, and it 
plays a primary role in many aspects of atmospheric 
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physics, including processes involved in meteorology, 
climate, the Earth’s radiation budget, and the global 
hydrological cycle. In addition, H 2 0 can be used as a 
tracer in studies of atmospheric transport and chem- 
istry. Therefore measurements of the spatial and 
temporal distribution of H 2 0 will lead to a better 
understanding of the various dynamic and chemical 
processes that occur in the atmosphere. An air- 
borne differential absorption lidar (DIAL) can make 
the necessary high-resolution measurements of H 2 0 
and aerosols over large regions of the lower atmo- 
sphere. This paper describes the design, develop- 
ment, and application of an airborne H 2 0 DIAL 
system for the remote detection of H 2 0 and aerosol 
profiles in the atmosphere. 

Several researchers have developed and imple- 
mented ground-based lidar systems to determine the 
vertical concentration profiles of H 2 0 remotely. 1-5 
These facilities include both Raman lidar, which 
incorporates inelastic backscattering to generate a 
return signal that is shifted in frequency from the 
frequency of the transmitted signal, and DIAL, in 
which elastic backscattering is combined with absorp- 
tion to produce a range-resolved profile. The capabili- 
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ties of many ground-based lidar systems have been 
enhanced by having them mounted on mobile plat- 
forms for measurements in diverse climates and 
under various atmospheric conditions. 6 8 To in- 
crease the flexibility of the lidar measurements fur- 
ther and to provide for mesoscale observations, one 
can design lidar systems to operate on airborne 
platforms. Because of the low-level return signals 
and the long averaging times associated with Raman 
lidar measurements, the method that is most adapt- 
able to performance of high-resolution measure- 
ments of H 2 0 from an aircraft is the DIAL technique. 5 
Thus all the current airborne lidar systems measure 
H 2 0 with the DIAL technique. 910 The NASA Lan- 
gley Research Center’s (LaRC) airborne H 2 0 DIAL 
system was designed to conduct mesoscale investiga- 
tions of H 2 0 and aerosol distributions under diverse 
meteorological conditions during both daytime and 
nighttime operation. 

Measurement Methodology 

A comprehensive approach was taken in the H 2 0 
DIAL research program to ensure that the airborne 
DIAL system would be capable of providing high- 
accuracy H 2 0 measurements over a large range of 
atmospheric conditions. A detailed error analysis of 
the DIAL technique for the measurement of H 2 0 
profiles was performed, and this analysis included 
random errors caused by lidar signal uncertainties 
and systematic errors caused by uncertainties in our 
knowledge of the H 2 0 absorption cross section. 
Experimental measurements were also conducted to 
improve our knowledge of the H 2 0 line parameters. 

H 2 0 Concentration Measurement 

DIAL Technique 

Because many authors have already presented thor- 
ough treatments of the DIAL technique, 11 12 only an 
overview of the theory is reviewed here. The DIAL 
procedure requires that one use high-energy, narrow- 
linewidth, tunable lasers to achieve accurate range- 
resolved gas measurements. One laser is tuned to 
an absorption line of the molecular species of interest 
(on-line), and the other is tuned to a nearby spectral 
region that has little or no absorption by the gas 
(off-line). The on-line and the off-line pulsed laser 
beams are transmitted through the atmosphere, and 
a portion of the laser radiation is backscattered by 
molecules and aerosols to a telescope receiver system 
that is colocated with the laser. The backscattered 
power from range R at the detector is given by 

PiR) = P Q { ct/2 )R 2 At\ (3 ( R ) 

x expj f | im{R) + k(R)\dR , (1) 

l*M) H 

where P 0 is the transmitted power, c is the velocity of 
light, t is the laser-pulse duration, R is the range, A is 
the receiver area, iq is the receiver-detector efficiency, 
P (R) is the atmospheric volume backscatter coeffi- 


cient, (j is the absorption cross section of the molecu- 
lar species of interest, n(R) is the number density 
profile of the gas, and k(R) is the atmospheric extinc- 
tion coefficient resulting from all other attenuation 
processes. 

Provided that the on-line and the off-line transmit- 
ters are close in wavelength (A A. < 0.1 nm), we can 
assume that p on {R) = p of f (/? ) and that k (m {R) = k { ^ (R), 
and we can solve Eq. (1) for the number density 
profile of the gas: 

1 P<tf{R'2}Pc>n(R\) o 

n[R) = 2(R 2 - R x )(<T (in - <r ofr ) ln P l ,^R l )P,jR 1 ) ' (2) 

where R 2 - R\ is the range cell for the average gas 
concentration, <r on - (r off is the differential absorption 
cross section for the two wavelengths, P on is the 
power received from range R for the on-line wave- 
length, and P ()fI is the power received from range R for 
the off-line wavelength. One can then convert the 
number density profile to a volume (or mass) mixing 
ratio by dividing the gas number (or mass) density by 
the ambient atmospheric number (or mass) density. 

DIAL Measurement Errors 

Errors in the measurement of H 2 0 with the DIAL 
method arise from both random and systematic 
sources. A discussion of these measurement errors 
for DIAL systems is available in the literature, 1,1-15 
and a detailed discussion of random and systematic 
errors for H 2 0 DIAL systems is given in Ref. 16. 
Random errors are caused by noise in the signal. 
Variability of the photon statistics in the lidar return 
signal, noise resulting from detector dark current, 
and noise in the background signal are the main 
sources of random error. One can minimize the 
random error associated with DIAL measurements 
by maximizing the signal-to-noise ratio (S/N). The 
S/N increase can be achieved in several ways, includ- 
ing an increase in laser-pulse energies or pulse repeti- 
tion rates, an improvement in the optical throughput 
of the system, the use of a detector system with 
higher quantum efficiency or lower dark-current noise, 
and the selection of a narrow-band optical filter that 
rejects most of the day background light and retains 
high optical efficiency. Following acquisition of the 
lidar data, we minimize random errors in the DIAL 
measurement by averaging the data. As a first step 
in noise reduction, the signals are averaged by smooth- 
ing them along the range, followed by averaging of a 
number of lidar profiles to enhance the precision of 
the measurement. This averaging of signals also 
results in a reduction in vertical and horizontal 
resolution. Thus a trade-off is necessary to achieve a 
balance between spatial resolution and measurement 
precision (typically a vertical range resolution of 
200-300 m and an averaging of 300-600 laser shot 
pairs are used in data reduction). Systematic errors 
in DIAL measurements are caused by both atmo- 
spheric and instrumental effects. The DIAL system 
presented here is designed and operated to minimize 
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the influences of systematic effects. Selection of the 
alexandrite on-line laser with a narrow linewidth and 
a high spectral purity, and its operation at the center 
of H 2 0 lines, ensures minimum influences in the 
DIAL measurement that are caused by the laser 
spectral distribution. The DIAL system is operated 
to avoid signal-overload effects, and the data analysis 
is restricted to regions that are not significantly 
influenced by signal overloads, signal nonlinearities, 
signal-digitization errors, and on-off signal-synchro- 
nization influences. The H 2 0 line selection ensures 
that negligible influence results from sensitivity of 
the H 2 0 cross section to uncertainity in the knowl- 
edge of atmospheric temperature. It is estimated 
that the sum of these systematic influences in cleaner 
regions of the lower troposphere will result in DIAL 
systematic errors of less than 5%. 

H 2 0 Spectroscopy 

Clearly, from Eq. (2), uncertainties in the knowledge 
of the gas absorption cross section can lead to errors 
in the DIAL measurements of gas profiles. Atmo- 
spheric influences on the effective cross section 
(i.e., the absorption cross section at the laser wave- 
length) are caused by the dependence of the absorp- 
tion linewidth and the absorption line-center position 
on pressure and temperature and by the dependence 
of the line strength on temperature. A high-resolu- 
tion spectroscopy experiment that incorporated a 
narrow-lindwidth tunable dye laser and two multi- 
pass absorption (White) cells with different pressure 
or temperature conditions has been used to provide 
measurements of H 2 0 absorption-line parameters in 
the 720-nm-wavelength region. 17 * 18 The use of the 
two multipass cells in parallel ensures that the results 
are insensitive to any laser-frequency drift during the 
measurement. Measurements of line strengths, pres- 
sure broadening, pressure-induced line shifts, and 
broadening and shift temperature dependences in the 
range 300-400 K were made. The pressure-shift 
coefficients in air were found to vary from line to line, 
and the average value was found to be -0.0170 
cm Vatm. The temperature exponent for the air 
broadening was found to be J dependent, with the 
broader lines (i.e., low-J lines) having the higher 
exponent. The average exponent value was mea- 
sured to be 0.67. The temperature-dependence expo- 
nent for the line shift in air was also found to be J 
dependent, but it showed the opposite behavior (i.e., 
the high-J lines had the higher exponents). The 
temperature exponent for the line shifting varied 
between 0.4 and 1.2. 

Using these laboratory measurements, we accu- 
rately determined the effective differential absorption 
cross section for atmospheric H 2 0, with an absolute 
error of 3—6%, for more than 120 H 2 0 absorption 
lines in the 720-nm region. These H 2 0 lines have a 
wide range of absorption cross sections and ground- 
state energies that can be selected to optimize the 
DIAL measurement according to the various atmo- 
spheric conditions encountered during the field experi- 
ments. 


Scattering-Ratio Measurement 

Along with the DIAL H 2 0 profile measurement, The 
DIAL system simultaneously provides information 
on the atmospheric aerosol distribution. An atmo- 
spheric scattering ratio is used to quantify the distri- 
bution of aerosols in the path of the laser beam. The 
scattering ratio is defined 11 as the ratio of the sum of 
the volume backscatter coefficients for the molecu- 
lar and the aerosol components of the return to the 
coefficient for the molecular signal only: 


R(r) = 


Pm(r) + Pfl(r) 

P 


1 + 


p» 

P Jr) ' 


(3) 


To determine the atmospheric scattering ratio, a 
measurement is first made in a clean region (0 Q ~ 0) 
to establish the relative molecular scattering profile. 
We then calculate the calibration constant by taking 
the ratio of the clean-region lidar return signal to the 
National Meteorological Center (NMC) molecular den- 
sity at that location. The calibration ratio is then 
used to normalize the NMC molecular density profile 
at any other location. This calculated molecular 
lidar return can be divided into the total lidar return 
to give the atmospheric scattering ratio. 


Airborne H 2 0 DIAL System 

In its present configuration, 9 the NASA airborne H 2 0 
DIAL system uses advanced instrumentation and 
technologies to provide enhanced H 2 0 measurement 
capabilities over the previous dye- laser-based H 2 0 
DIAL system. 19 20 The system has evolved over the 
past 5 years to include tunable solid-state laser 
technology; equipment for accurate wavelength posi- 
tioning, stabilization, and monitoring; advanced data- 
acquisition and experiment-control electronics; high- 
throughput receiver optics; and solid-state detector 
technology. The airborne H 2 0 DIAL system con- 
sists of five major subsystems: the laser transmit- 

ters and transmission optics, the equipment for wave- 
length positioning and monitoring the receiver optics 
and detectors, the timing and control electronics, and 
the data-acquisition instrumentation. The laser 
transmitters and the transmission optics are con- 
tained, along with the wavelength-monitoring equip- 
ment, in a rigid structure that is bolted to the aircraft 
seat tracks. The receiver system is also attached to 
this structure, and it can be mounted for either nadir 
or zenith operation. Three dual aircraft racks house 
the data-acquisition and the experiment-control elec- 
tronics, the laser power supplies and heat exchangers, 
and the wavelength-control electronics. A fourth 
rack, smaller than the dual racks, holds the two color 
printers used for real-time data output. A schematic 
of the airborne H 2 0 DIAL system in the NASA 
Wallops Flight Facility (WFF) Electra aircraft is 
shown in Fig. 1. 

Laser Transmitters 

Different laser systems are used to generate the 
on-line and the off-line wavelengths for the H 2 0 
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Fig. 1 . Schematic diagram of the NASA LaRC airborne H 2 O DIAL system. 


DIAL measurement. A tunable, narrow-band alex- 
andrite laser is used as the on-line transmitter to 
provide the requisite linewidth, wavelength stability, 
and spectral purity for accurate measurement of H 2 0 
concentration. Because the spectral requirements 
for the off-line laser are less stringent than those for 
the on-line laser, an Nd:YAG laser-pumped dye laser 
is used to generate the unabsorbed wavelength. 

On-Line Laser 

The alexandrite laser used as the on-line transmitter 
(Fig. 2) is a custom-designed laser system developed 
by Allied Military Laser Products. A double-ellipse 
pump chamber contains two flashlamps and a 110 
mm x 5 mm (length x diameter) alexandrite rod with 
antireflection-coated plano-plano ends. Because the 
alexandrite emission cross section increases with 
temperature, the laser rod is heated to 70 °C with a 
heat exchanger to optimize the laser gain. A sepa- 
rate cooling loop is used to operate the lamps at a 
lower temperature of 40 °C. The Fabry-Perot cavity 


that encloses the pump chamber is composed of a 
reflector with a 5-m-high radius of curvature and a 
60%-reflectivity piano output coupler. This resona- 
tor cavity is Q switched with an acousto-optic device 
to produce a 200-ns pulse length at a rate of 10 Hz. 
The laser wavelength and linewidth are controlled 
with a birefringent tuner for coarse control and a pair 
of Fabry-Perot etalons for fine control. The birefrin- 
gent tuner, which consists of a stack of five quartz 
plates in a ratio of 1:2:2:10:10, is angle tuned with a 
stepper motor. A 1-mm-thick plate of fused silica 
with 30%-reflectivity coatings is used for the low- 
resolution etalon. This element serves as an order 
selector for the high-resolution etalon and is also 
angle tuned by a stepper motor. The high-resolu- 
tion etalon is a 1-cm air-spaced Fabry-Perot etalon 
that consists of two wedged, fused silica plates with 
30%-reflectivity coatings on the inside surfaces and 
antireflection coatings on the external surfaces. We 
tune this element by varying the spacing of the plates 
with piezoelectric ceramics. Synchronous scanning 
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Fig. 2. Optical configuration of the alexandrite laser cavity. 
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of the laser over a sufficient wavelength interval 
without mode hopping requires that the free spectral 
range ratio between the two etalons be very close to a 
half-integer value. This ratio is adjusted by use of a 
scanning single-longitudinal-mode ring dye laser as a 
calibration source. 21 The three intracavity tuning 
elements together provide an output linewidth of less 
than 1.1 pm (FWHM). They can be synchronously 
scanned through a microprocessor-controlled routine 
to vary the laser wavelength continuously over a 
150-pm region located anywhere in the 725-785-nm 
tuning range of the laser. Because the large number 
of coated intracavity elements, the output pulse 
energy of the laser is limited to 30 mJ to minimize the 
possibility of optical damage. 

For accurate DIAL measurements, the laser line 
must be nearly monochromatic with respect to the 
absorption line; i.e., it must have high spectral purity 
and a narrow linewidth. Spectral purity is defined 
as the percentage of laser energy transmitted within a 
specific spectral interval relative to the total transmit- 
ted laser energy. For the H 2 0 lines in the 730-nm 
region, the spectral purity must be > 99% to achieve 
high-accuracy H 2 0 measurements. 16 In addition to 
high spectral purity, the laser linewidth must be 
small compared with the linewidth of the H 2 0 absorp- 
tion feature. Finite laser linewidths lead to uncer- 
tainties in the knowledge of the absorption cross 
section and cause systematic errors in the DIAL 
measurements of the H 2 0 profiles. 1516 

The spectral properties of the alexandrite laser 
were measured with a multipass absorption cell and a 
ratiometric energy meter. 22 ’ 23 These multipass-cell 
transmission values indicated that the spectral purity 
of the alexandrite laser was >99.85%. To deter- 
mine the influence of the finite laser linewidth, 
measurements of the effective absorption cross sec- 
tion were made for different values of the absorption 
linewidths of H 2 0 and oxygen. For absorption line- 
widths greater than 4 pm (FWHM), which correspond 
to altitudes below 10 km, the cross section is within 
1% of the cross section calculated from spectroscopic 
data. 18 22 * 23 Therefore, under these conditions, the 
alexandrite laser can be considered to be monochro- 
matic with respect to the H 2 0 absorption line. 

In addition to having a monochromatic spectral 
distribution, if the line-center absorption cross sec- 
tion is used in the data-reduction process, a precise 
positioning of the laser wavelength relative to the 
H 2 0 absorption line center is required for accurate 
measurement of H 2 0 concentration. To obtain a 
DIAL H 2 0 measurement error of less than 5% at 
altitudes below 5 km, the laser-wavelength deviation 
from the center of the absorption line must be less 
than ±0.8 pm, including initial positioning uncer- 
tainty, shot-to-shot wavelength jitter, and long-term 
wavelength drift. 16 Without additional wavelength 
stabilization of the alexandrite laser, however, mode 
hopping and wavelength drifts as large as 3 pm h 
have been observed in a controlled laboratory environ- 
ment. Consequently, a stabilization device was de- 


signed and implemented to maintain the required 
spectral characteristics despite the pressure, humid- 
ity, and temperature variations and the mechanical 
vibrations inherent in an aircraft environment. 22 23 

A Laboratory for Science model 200 wavelength- 
stabilized cw He-Ne laser is used as an external 
absolute reference to stabilize the optical thickness of 
the two Fabry-Perot etalons actively. We generate 
two error signals by comparing the output intensity 
of the He-Ne laser with the intensities of the He-Ne 
laser beam that reflects from the air-gap etalon and 
that transmits through the solid etalon. The angle 
of the He-Ne beam through each etalon is adjusted so 
that the He-Ne laser wavelength corresponds to the 
half-maximum of each of the Fabry-Perot transmis- 
sion modes. At this point there is a maximum 
sensitivity to changes in the optical thickness of the 
etalons. The three intensity levels are digitized and 
sampled by a computer, and the computer creates the 
appropriate correction voltages, which are applied to 
the etalon controllers. Using an absorption feature 
of oxygen contained in a multipass cell as a reference, 
a long-term ( ~ 1.5 h) laser drift of less than 0.7 pm 
relative to the absorption-line center was observed. 

Off-Line Laser 

As shown in Fig. 3, the off-line laser wavelength is 
generated by use of a Quantel model 482 frequency- 
doubled Nd:YAG laser to pump a Jobin-Yvon model 
HP-HR tunable dye laser. 19 20,24 The Nd:YAG laser 
oscillator uses a stable resonator configuration to 
produce a TEM 0 o transverse-mode output beam at 
1.06 jxm, with a pulse length of 12 ns. After being 
up-collimated by a telescope, the oscillator beam is 
passed through two amplifiers, giving it an output 
energy of 1 J/pulse at a repetition rate of 10 Hz. For 
frequency doubling to 532 nm, an angle-tuned, tem- 
perature-stabilized KDP-II crystal with 40% effi- 
ciency is used to generate 400 mJ/ pulse. Two di- 
chroic beam splitters are used to separate the visible 
and the infrared wavelengths and to direct the green 
pump beam into the dye laser. 

To produce the high-energy, narrow-linewidth out- 
put pulses required for accurate DIAL measure- 
ments, the dye laser also uses an oscillator-ampli- 
fier configuration. Narrow-linewidth operation is 
achieved with incorporation of a Littrow-mounted 
holographic grating in the oscillator cavity. In com- 
bination with an intracavity Pellin-Broca prism, the 
grating (2400 lines/ mm) produces a linewidth of less 
than 15 pm, which is adequate to satisfy the off-line 
laser requirements. The line-narrowed output from 
the oscillator is amplified in three stages: the first 

stage is a preamplifier that operates below saturation, 
and the final two stages are amplifiers that operate 
near saturation. By sufficient saturation of the am- 
plifiers, the input energy is efficiently transferred to 
the narrow-band portion of the laser beam, and 
amplified spontaneous emission is minimized. 20 
The amplified spontaneous emission generated by the 
preamplifier is suppressed by two Pellin-Broca prisms 
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placed between the preamplifier and the first amplifier. 
The dye solution used for operation near 730 nm is 
Exciton LDS-751 in methanol. With this dye solu- 
tion the laser produces an output energy of 60 mJ, 
which corresponds to a conversion efficiency of 15%. 
To permit the use of the same detection system in the 
receiver for the on-line and the off-line lidar returns, 
the two laser systems are fired with 300-gs separation. 
This is close enough to sample the same atmosphere 
with the laser beams but far enough apart to have 
independent lidar return signals. In addition, even 
though the 200-ns pulse length of the flashlamp- 
pumped alexandrite laser is much longer than the 
12-ns pulse of the dye laser, this can be compensated 
for in the data-processing procedure. Performance 
characteristics of both the on-line and the off-line 
airborne H 2 0 DIAL laser transmitters are summa- 
rized in Table 1, and a block diagram of the laser 
transmitters, wavelength-control equipment, and 
spectral-monitoring equipment is presented in Fig. 4. 

Transmission Optics 

Before they are transmitted into the atmosphere, the 
alexandrite and the dye-laser output beams pass 
through various optical components for attenuation, 
evaluation, and steering. First, as a result of ground 
observer eye safety considerations during nadir opera- 
tion, each beam is sent through a variable attenuator 
consisting of two polarizers with a A/2 plate between 
them. Depending on aircraft altitude, we adjust the 
output beam energies with the attenuators to achieve 
a combined radiant exposure at ground level that is 
below the maximum permissible exposure value as 
determined by the American National Standards 
Institute. 26 Energy monitors, which use less than 
2% of the output laser energy, are mounted after the 
attenuators, and the shot-to-shot output energy is 
measured and recorded. The beams are then re- 
flected off of separate mirrors held in motorized 
gimbal mounts for independent optimization of each 
beam in the receiver field of view. Both beams are 
finally reflected off a single mirror mounted behind 


the telescope secondary and transmitted coaxially 
with the receiver telescope into the atmosphere. 
The beams are completely enclosed in baffles until 
they reach the aircraft window to prevent scattered 
light from entering the telescope. Beam steering 
and energy attenuation are accomplished remotely 
from the data-acquisition system. 

Monitoring Instrumentation 

A McPherson model 2051 spectrometer and a Spectra 
Physics model LO-3 absorption cell (see Fig. 4) are 
used to position the on-line transmitter wavelength 
at the center of an H 2 0 line and to monitor the 
performance of the laser during flight. The 1-m 
spectrometer incorporates a 1200 lines mm holo- 


Table 1 . Airborne H 2 0 DIAL Transmitter Characteristics 


Characteristic 

Value 

On-Line 

Alexandrite laser 

(Allied model 407000) 

Wavelength 

72.5-780 nm 

Pulse energy 

30 mJ 

Pulse length 

200 ns 

Repetition rate 

10 Hz 

Linewidth 

1 . 1 pm 

Spectral purity 

> 99.853 

Wavelength stability 

0.7 pm 

Wavelength jitter 

< 1 0.5 pm 

Off-Line 

Nd:YAG laser 

jQuantel model YG482I 

Wavelength 

532 nm 

Pulse energy 

400 mJ 

Pulse length 

15 ns 

Repetition rate 

10 Hz 

Dye laser 

Uobin-Yvon model IIP- HR) 

Wavelength 

726-730 nm 

Pulse energy 

30 mJ 

Linewidth 

' 15 pm 
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Fig. 4. Block diagram of the H 2 0 DIAL laser transmitter sub- 
system including laser transmitters, wavelength-control equip- 
ment, and spectral-monitoring instruments. JY, Jobin-Young, 
OMA, optical multichannel analyzer. 


graphic grating in a Czerny-Turner optical configura- 
tion to produce a linear dispersion of approximately 
0.7 nm/mm at 730 nm. A Princeton Applied Re- 
search optical multichannel analyzer (OMA), which 
consists of a model 1205D linear array detector head 
attached to the spectrometer and a model 1205 A 
control unit, is used to evaluate the spectral character- 
istics of the output beam. The detector has a pixel 
spacing of 25 pm, which corresponds to a wavelength 
dispersion of 1 7.5 pm / pixel and provides a wavelength- 
positioning accuracy of better than ±20 pm. 

The multipass absorption cell is of stainless-steel 
construction and has a heating jacket for high- 
temperature operation. The cell mirrors have dielec- 
tric coatings with > 99% reflectivity over the range 
725-790 nm. During operation, the cell is evacu- 
ated, heated to 60 °C, and filled with approximately 40 
Torr of pure H 2 0. Because the saturation vapor 
pressure at 60 °C is 150 Torr, the risk of condensation 
on the mirrors is minimized. A Laser Precision 
model RJ7200 ratiometric energy meter with model 
RJP 765 silicon probes is used to measure the amount 
of absorption through the multipass cell. A stability 
test was performed to determine the ratiometer 
precision associated with this type of measurement. 
The ratiometer precision, defined as the variation in 
the absorption-cell transmission divided by the aver- 
age absorption-cell transmission, was determined to 
be 1 % when we evaluated the reproducibility of the 
beam transmission through the evacuated absorption 
cell with the laser set at a fixed wavelength. 23 
Taking into account the ratiometer precision, there is 
sufficient optical depth, with an H 2 0 pressure of 40 
Torr and a path length of 100 m in the cell, to position 
the alexandrite wavelength accurately at the center of 
H 2 0 lines having absorption cross sections greater 
than 4 x 10 -24 cmVmolecule. 

Thus the procedure implemented to optimize the 
alexandrite output on an H 2 0 line is first to set the 
wavelength as close as possible to the selected H 2 0 
line by use of the spectrometer-OMA, with an atomic 


vapor lamp as a reference. Next, while observing the 
absorption-cell transmission on the ratiometer read- 
out, we vary the laser wavelength by synchronously 
scanning the intracavity tuning elements. When the 
maximum absorption through the multipass cell is 
detected, the laser wavelength is locked on the H 2 0 
line by use of the stabilization device described above. 
Using the spectrometer-OMA combination, we set 
the off-line wavelength between 60 and 100 pm from 
the on-line wavelength in a region free of other 
absorption features. 

Both the spectrometer-OMA combination and the 
absorption cell are also used for real-time monitoring 
of the alexandrite spectral characteristics during 
flight operations. A small portion of the output is 
kicked off and sent through fiber optics to the spec- 
trometer. Mode hops between adjacent thin etalon 
modes (free spectral range 180 pm) that result from a 
drift of the birefringent tuner or of the etalon itself 
can be seen on the OMA display. A beam splitter is 
also placed in the beam to reflect 4% of the output 
into the absorption cell. Thus any drift in wave- 
length that is caused by a variation in the optical path 
through the air-spaced etalon will be immediately 
detected by the ratiometer. In addition, a shot-to- 
shot measurement of the effective absorption cross 
section in the absorption cell can be recorded. 

Receiver System 

Figure 5 shows the present optical configuration of 
the receiver system. The receiver system incorpo- 
rates a 36-cm-diameter Celestron model C-14 tele- 
scope with Schmidt-Cassegrain optics to collect the 
backscattered laser light and focus it into the detector 
package. An iris placed in the telescope focal plane 
acts as a variable field stop for adjusting the field of 
view of the receiver. The field of view is normally set 
at 2 mrad to contain greater than 97% of the energy 
in the approximately 1.5-mrad-divergence laser beam. 
In addition to having interchangeable detectors, the 
detector package contains optics for beam collima- 
tion, filtering, and focusing. After passing through 
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Fig. 5. Optical configuration of the H 2 0 DIAL receiver system. 


6428 


APPLIED OPTICS / Vol. 33, No. 27 / 20 September 1994 










the field stop, the laser light is collimated by a 
100-mm-focal-length lens and prefiltered with a di- 
chroic beam splitter. The lidar signals are then 
transmitted through an interference filter to reject 
background light and any Raman scattered laser light 
further. 27 After it is transmitted through the filter, 
the laser light can be directed to two possible detector 
positions. In addition, the available return power 
can be split between a near-field detector with low 
gain and a far-field detector with high gain to improve 
the dynamic range of the measurement. 24 Various 
filter and detector configurations were used in flight 
tests conducted during the field measurement cam- 
paigns, and substantial improvements were imple- 
mented between each field deployment. Table 2 lists 
the important characteristics of the airborne H 2 O 
DIAL receiver system. 

Receiver Filters 

During daytime operation, a 0.4-nm-bandwidth 
(FWHM) filter is required for high solar-background 
rejection. The daytime filter (Spectro-Film 7289/3-4 
SF) used during the 1989 and 1990 field measure- 
ments had a relatively low transmission of 33%. 
However, a new filter (Andover ANDV1551) with a 
transmission of 48% was acquired before the 1991 
flight tests. This represents a receiver system effi- 
ciency increase of 45% and a resulting measurement 
S/N increase of 20%. The increase in S/N is re- 
flected in the 1991 H 2 0 data, which required substan- 
tially less signal shot averaging than the 1989-1990 
data. At night, a higher-transmission 1-nm-band- 
width filter (Spectro-Film 7289/10 SF) replaces the 
daytime filter, beause only the Raman-backscattered 
component needs to be blocked. In addition, night- 
time measurements can be made with no filter in 
place if corrections are made for the Raman return 
signal, which is approximately 3% of the Rayleigh 
return signal. 27 The filters are placed in an oven and 
heated to 40 °C to stabilize their transmission peaks 
thermally at the on-line wavelength. The H 2 0 lines 
in a 4-nm-wavelength range can be achieved with 
each filter by use of angle tuning with a motorized 
encoder rotation stage. 

Receiver Detectors 

During the field deployments, a photomultiplier tube 
(PMT) and an avalanche photodiode (APD) were 
incorporated as detectors for both daytime and night- 


Table 2. Airborne H 2 0 DIAL Receiver Characteristics 

Characteristic Value 


Receiver area 
Field of view 
Filter bandwidth 
Filter transmission 
Optical efficiency 
Detector type 
Quantum efficiency 
Total efficiency 


0.1 m 2 
1.5 mrad 
0.4 nm 
48% 

34% 


PMT 

APD 

3.9% 

85% 

1.3% 

29% 


time measurements in diverse meteorological condi- 
tions to determine which detector exhibited optimum 
performance in the system. The choice of a detector 
is very important because it can influence both sys- 
tem throughput and dark current. An APD detector 
provides high quantum efficiency, but it has increased 
noise because of high dark current and a large 
excess-noise factor. 16 Compared with an APD, a 
PMT has lower quantum efficiency, but at the same 
time the PMT has a lower dark current and a lower 
excess-noise factor. As has been discussed in the 
literature, 16 28 29 the APD is more suitable for DIAL 
measurements in cases in which signals are strong 
compared with the detector dark current, and a PMT 
is more suitable for measurements in low-signal 
regimes. 

The RCA model 7265 PMT has an enhanced S-20 
photocathode for improved operation in the red and 
the near-infrared regions of the spectrum. The AP D 
module includes a 1.5-mm-diameter RCA model 
C30955E APD and its preamplifier, motorized en- 
coder translation stages to optimize the position of 
the APD in the field of view, and focusing optics to 
image the field stop on the APD surface. The APD 
data indicated a much higher measurement S/N, 
thus leading to higher precision and longer-range 
measurements with this detector. This conclusion 
is evident from Fig. 6, which shows a comparison of 
the measured H 2 0 mass mixing ratios (M ), the mix- 
ing ratio standard deviation (<r), and the measure- 
ment signal-to-noise ratio (M/a) from a nighttime 
flight test on 18 April, 1991. We accomplished the 
detector comparison by conducting measurements on 
a specific flight path with the PMT and, after chang- 
ing detectors during the turn, retracing the same 
path for measurements with the APD. The data 
shown in Fig. 6 have been calculated with a 210-m 
vertical range cell size and averaged in the horizontal 
direction over approximately 12 km. The S N plots 
clearly indicate the superior performance of the APD 
over the entire altitude range, with the exception of 
the final 300 m of the measurement. Improvement 
over the PMT is expected to be even greater during 





(c) SIGNAL TO NOISE RATIO 


(b) STANDARD DEVIATION 

Fig. 6. Performance comparison of the APD and the PMT 
detectors from a flight on 18 April 1991. The plots show (a.i the 
measured mixing ratio profiles, (b) the mixing ratio standard 
deviation profiles, and (c) the mixing ratio S/N profiles. AGL, 
above ground level. 
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daytime operation because the PMT S/N is degraded 
by background radiation, whereas the APD S/N is 
still dark current limited. 29 

Timing, Control, and Signal-Processing Electronics 

As mentioned above, the on-line and the off-line laser 
pulses are transmitted sequentially, with a temporal 
separation of 300 gs. This approach reduces re- 
ceiver system complexity by eliminating the need to 
separate the DIAL wavelengths spectrally. How- 
ever, the simplification of the receiver system is at the 
expense of increased complexity in the design of the 
timing and control electronics. 

The inherent characteristics of the acousto-optic Q 
switch used in the alexandrite laser further compli- 
cate the requirements of the timing and the control 
functions. For example, the output from the alexan- 
drite occurs approximately 4.7 ps after the lase 
command is issued and has a 150-ns jitter, whereas 
the output pulse from the Nd:YAG system occurs 
0.23 gs after its lase command, with a jitter of only 
3 ns. 

The most important factor in the timing methodol- 
ogy is the synchronization of the sequential laser 
firing with the return signal digitization. To fulfill 
the timing requirements, a master control unit was 
developed to control laser flashlamp firing, lase fire 
command, and PMT gating with a precision of ± 10 
ns. In addition, a lase-coherent time base was de- 
signed to provide a digitizer trigger that is synchro- 
nized with the actual laser firing, which is detected 
with fiber-optic pickup by use of a fast photodiode 
incorporated in the time base. The time base of the 
digitizer is stopped and restarted with the correct 
phase relative to the time of the laser firing for both 
the on-line and the off-line laser pulses. This phase 
correction was found to be necessary because an 
improper delay between the on and the off signals can 
distort the DIAL-measured profile in regions of strong 
aerosol gradients. 30 

The lase-coherent trigger to the digitizer is gener- 
ated at the time of the alexandrite laser firing, and 
any timing jitter between the alexandrite and the 
Nd:YAG laser pulses would result in an apparent 
desynchronization of the dye-laser return in the 
digitizer memory. Because the alexandrite laser does 
have a significant amount of jitter associated with it, 
a modification was implemented to create a precision 
delay in the Nd:YAG output referenced to the on-line 
laser pulse by use of a Data Design model DG11A 
digital delay generator. Other features include 
phase-aligned, precision-delayed lase markers summed 
with the return signals and line-locked firing-rate 
logic that minimizes 60-Hz drift in the dc offset. 

When a PMT is used as the detector in the receiver 
system, it is gated on by a transistor-transistor logic 
pulse with an adjustable time delay after each laser 
pulse is transmitted. PMT gain is controlled by 
manual ajdustment of the high voltage supplied to 
the PMT dynode network. The APD, however, is a 
nongated detector with fixed gain. The return sig- 


nals are filtered (2.5-MHz bandwidth) and amplified 
by a factor of 10 with a Digital Signal Processing 
(DSP) model 1020 amplifier. The signals are then 
digitized at 10 MHz by use of a 12-bit DSP model 
2112S transient recorder with averaging memory. 
Both the on-line and the off-line DIAL returns are 
processed sequentially during a 4096-Word sweep of 
the digitizer. Each sweep can then be stored and 
summed in the averaging memory. Typically, 15 
sweeps are summed before they are transferred to the 
host computer. The digitizer and the amplifier are 
interfaced to a computer automated measurement 
and control (CAMAC) crate, and certain functions, 
such as the number of sweeps summed, are directly 
programmable from the computer. 

Data-Acquisition System 

The DIAL Data Acquisition System (DAS) is an 
advanced version of the system described in Ref. 31. 
The DAS occupies a standard aircraft double rack 
with digitizers, control electronics, and PMT power 
supplies in one half and the computers and peripher- 
als in the other half. The digitizer and the amplifier 
are mounted in a CAMAC crate, and they communi- 
cate with the rest of the data system through a 
Kinetics System model 3922 crate controller. 
Through this crate controller the computer system is 
able to set up various options as well as initialize and 
enable the digitizer. When the programmed number 
of sweeps has been accomplished, an interrupt is 
generated by the digitizer, and the data are then 
transferred from the averaging memory to the host 
computer through the crate controller. Other in- 
puts to the computer system include laser energies, 
aircraft position data ( latitude and longitude informa- 
tion from the aircraft inertial navigation system), 
pressure, total temperature, dew-point temperatures, 
and absorption-cell ratiometer readings. Laser ener- 
gies are detected and converted from analog to digital 
through a custom-designed energy monitor system. 
The aircraft position data are in an ARINC-561 
format and are decoded by custom-designed electron- 
ics. Energies and navigation data are transferred to 
the host computer through two 16-bit parallel inter- 
faces. Pressure, total temperature, and dew point 
from two instruments are converted from analog to 
digital signals, sent through light pipes from the front 
of the aircraft to the DIAL system at the rear, and 
transferred to the computer through a third parallel 
interface. The ratiometer reading of the absorption- 
cell transmission is converted to RS232 format and 
interfaced to the computer through a four-channel 
serial input-output CAMAC module. 

The DIAL DAS is based on two Digital Equipment 
Corporation (DEC) PDP 11/73 processors and one 
MDB JFEP-11 front-end processor. Each PDP has 
128 kWords of 16-bit memory, although only 28 
kWords of memory is available for resident program- 
ming. The overall flow chart for the DIAL DAS is 
shown in Fig. 7. In general, all data acquisition and 
storage are performed by the PDP at the left in Fig. 7 
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Fig. 7. Block diagram of the H^O DIAL Data Acquisition System, 
operating system. 


INS, inertial navigation system; SLU, serial line unit; MOS, master 


(PDP-l), whereas all data display and analysis are 
performed by the PDP at the right in Fig. 7 (PDP-2). 
The JFEP-1 1 is housed in PDP-2 and has 256 kWords 
of memory, of which 128 kWords can be shared with 
the host computer. The DEC RT11-SJ operating 
system is used for software development. The oper- 
ating system and the DIAL software are stored on and 
retrieved from 300-Mbyte hard disks. Also available 
on each PDP is a 5-1/4-in. floppy disk drive for hard 
disk backup. These floppy disks are also used for 
data input to the DIAL DAS from external comput- 
ers. 

PDP-l acts as the master computer through which 
the operator communicates with PDP-2 and the 
JFEP. The operator communicates with the proces- 
sors through a Matrox CTM-300 smart keyboard, and 
input— output is echoed on a small black and white 
monitor. The keyboard interfaces to PDP- 1 through 
the console port of a four-port serial line interface 
(RS232). A second port of this interface connects to 
PDP-2, and a third port connects to the JFEP serial 
line unit. 

The DIAL DAS operating system software consists 
of three programs that run concurrently on the two 
PDP processors and on the JFEP processor. When 
the DIAL master program is being loaded on PDP-l, 
commands are automatically generated to load and 
run appropriate programs on the other two processors. 
In the data-acquisition mode, PDP-l will collect the 
digitizer data, the supplemental information men- 
tioned above, and the time and date from a crystal 
clock into one long data buffer. This buffer can 
represent anywhere from 1 to 65536 laser pairs 
summed in the digitizer averaging memory. Once 
the digitizers achieve the selected number of sweeps. 


an interrupt is sent to alert PDP-l. The auxilliary 
information is then attached in the form of a buffer 
header, and the entire buffer is then recorded on 
8-mm cartridge tape in binary format. The data 
buffer and a table of display options are then sent to 
PDP-2 through a direct-memory-access interproces- 
sor link. The data buffer and the display table are 
shared between PDP-2 and the JFEP through dual- 
ported memory. 

The user controls data-acquisition and display op- 
tions by keying in commands specifically designed for 
this DIAL system. The command structure format 
provides a fast update of user-selected options. 
Function keys are defined for some of the more 
frequently used commands. 

As stated above, PDP-l is mainly dedicated to data 
acquisition and storage and to communication with 
the operator. Additionally, if a particular display 
option is selected, PDP-l will process the data for 
aerosol information and compute either the relative 
aerosol scattering ratio or the total aerosol scattering 
ratio. The aerosol profiles are then color coded and 
presented as a function of time and altitude on one of 
the ACT II inkjet printers. 

The main function of the second PDP is to perform 
selected analysis options for display on a black and 
white monitor through a graphics display controller. 
The display modes are designed to provide diagnostic 
tools, feedback on beam optimization, and results of 
calculations of aerosol and H 2 O profiles. The most 
basic mode of display shows the returns from both 
lasers in an overlapped manner as a function of range. 
This display is used to observe signal levels during 
beam steering and beam optimization. The operator 
may elect to display only one of the returns, subtract 
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the background, smooth the return (by use of a 
running average), correct for range-squared depen- 
dence, normalize the return for laser energy, and 
average successive returns. Another display mode 
shows the relative aerosol scattering profile as a 
function of altitude with smoothing and averaging 
options. One can also calculate scattering ratios in 
this mode by specifying an altitude at which to 
normalize the lidar return to a model molecular 
number density profile. A H 2 0 mixing ratio profile 
can be calculated from each DIAL signal pair and 
successive profiles can be averaged to improve the 
measurement statistics. The updating average can 
be displayed as a function of altitude with range-cell 
size, smoothing, and averaging options. 

A final display mode is implemented during spiral 
ascents and descents. The H 2 0 mixing ratios are 
calculated from in situ dew-point measurements and 
plotted as a function of altitude. These profiles can 
then be compared with inflight DIAL profiles ob- 
tained over the spiral point to provide validation of 
the DIAL system. Any screen display can be cap- 
tured at any time and plotted with a dot matrix 
printer. Because the screen is stored in memory 
before plotting, there is only a brief pause before the 
screen update proceeds. 

PDP-2 sends the DIAL returns to the JFEP. The 
JFEP stores a running average of DIAL returns and 
computes the DIAL calculation on that average to 
obtain a H 2 0 mixing ratio profile. Each profile is 
then returned to PDP-2, color coded for magnitude, 
and plotted on an ACT II inkjet printer as a function 
of altitude, time, and position. Real-time color cross 
sections of H 2 0 mixing ratios and aerosol scattering 
ratios are produced simultaneously by the H 2 0 DIAL 
DAS. 

Field Measurements 

During 1989-1991, the H 2 0 DIAL system was de- 
ployed to NASA WFF for three separate flight-test 
experiments. Preliminary engineering flights were 
performed in July 1989 to complete system shake- 
down and to assess system airborne operational reli- 
ability. 32 During March 1990 and April 1991, addi- 
tional flight tests were conducted in a variety of 
meteorological conditions to verify the performance 
of the system and to establish the accuracy of the 
DIAL-measured H 2 0 profiles. 9 * 33 Included in these 
engineering and data flights were numerous intercom- 
parison measurements taken with onboard dew-point 
hygrometers and H 2 0 radiosondes launched from 
both NASA LaRC and NASA WFF. Between each of 
the field deployments, many modifications were imple- 
mented to improve system performance and to auto- 
mate system operation in flight further. 

In Situ Instruments 

Two dew-point hygrometers were used on the NASA 
WFF Electra aircraft to obtain in situ H 2 0 measure- 
ments: the EG&G model 137-C3 aircraft hygrom- 
eter and the modified EG&G model 300 microproces- 


sor-controlled hygrometer. Both instruments use 
the chilled-mirror dew-point condensation principle 
to determine the H 2 0 concentration outside the 
aircraft. Details of this technique and of the opera- 
tion of an EG&G 137-C3 hygrometer in aircraft field 
measurements can be found in Ref. 34. This unit 
incorporates a model S10 aircraft-mounted sensor 
head with a sampling probe designed to provide a 
static pressure sample at the correct mass flow rate 
across the chilled mirror. The model 300 hygrom- 
eter, on the other hand, is a laboratory unit, and we 
had to modify its model 53 sensor head by implement- 
ing a sampling probe for airborne operation. A 
major advantage of this instrument is the application 
of microprocessors and digital techniques to control 
the servoloop that maintains dew or frost in equilib- 
rium on the mirror surface in the sensor. This 
approach leads to greater accuracy as well as the 
ability to operate effectively at lower dew points. 
The accuracy of the model 300 hygrometer is specified 
as ±0.2 °C, whereas the accuracy of the model 137-C3 
hygrometer is stated to be ±0.5 °C. Conversion of 
the vapor pressures, determined with the dew-point 
temperature measurements, to H 2 0 mixing ratios is 
accomplished by making simultaneous measure- 
ments of the atmospheric pressure. In situ H 2 0 
mixing ratio profiles were obtained with the dew- 
point hygrometers by spiraling the Electra down to a 
minimum altitude of 152.4 m at approximately 304.8 
m/min. 

Other in situ intercomparison profiles were ac- 
quired by H 2 0 radiosondes launched from ground 
stations located at NASA LaRC and NASA WFF. 
The balloons were released approximately 10 min 
before overflight by ground personnel who were in 
radio contact. Because the balloons typically ascend 
at approximately 200 m/min, this timing places the 
sonde at the center of the DIAL measurement alti- 
tude range when the plane is overhead (neglecting 
horizontal movement resulting from winds). The 
radiosonde package used for these measurements was 
the Vaisala model RS80, which incorporates solid- 
state thin-film capacitor sensing technology. This 
type of sonde was chosen because of its sensitivity and 
accuracy at relative humidity values below 20% as 
compared with the carbon element hygristor tech- 
nique used in other packages. 35 Numerous over- 
flights of NASA WFF and NASA LaRC were conducted 
to provide ample opportunity for DIAL-radiosonde 
comparisons. In addition, a standard landing proce- 
dure was used at NASA WFF, first to overfly the 
facility to obtain a DIAL profile and then to spiral the 
aircraft down while simultaneously launching a H 2 0 
radiosonde. This landing procedure allowed us to 
intercompare all three H 2 0 profile measurement 
techniques. 

System Parameters 

Two H 2 0 absorption lines near 730 nm were used for 
the majority of the DIAL measurements. The stron- 
ger line, located at 728.7379 nm, has an absorption 
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cross section of 46.0 x 10 24 cm 2 /molecule at ground 
level. This line was used extensively in the flight 
tests because of the ability to make high-accuracy 
measurements from 4 km to the surface over a broad 
range of meteorological conditions. The second line, 
with an absorption cross section of 21.6 x 10 24 
cm 2 /molecule at 1 atm, provides approximately a 
factor-of-2 decrease in the optical depth over the 
range cell compared with the stronger line. Measure- 
ments from higher altitudes and at elevated humidity 
levels required that we use this weaker line, which is 
located at 728.8126 nm. These two H 2 0 absorption 
lines were chosen for three primary reasons: ( 1 ) The 

receiver daytime filter center wavelength for light at 
normal incidence is 729.2 nm; (2) errors resulting 
from temperature sensitivity are small because both 
lines have lower energy levels in the temperature- 
insensitive range for mixing ratio measurements 36 ; 
and (3) the lines are isolated, thus minimizing the 
errors caused by overlapping from adjacent lines and 
facilitating the positioning of the off-line wavelength. 
The dye-laser wavelength was placed at approxi- 
mately 728.90 nm when the weaker line was used and 
at approximately 728.65 nm when the stronger line 
was used. In addition to providing the off-line laser 
return for the DIAL measurement, the dye-laser 
backscattered signal is also used to produce the 
atmospheric scattering-ratio profile. 

Data Analysis 

The postmission data analysis was done with a Mi- 
crovax II computer. The data were first edited by 
rejection of shots resulting from poor laser energy or 
signals that were attenuated by strong near-field 
clouds to ensure that quality lidar return signals were 
used in DIAL calculations. The on-line and the 
off-line signal profiles were altitude aligned by shift- 
ing them to compensate for any slight difference in 
the altitude registering caused by a relative laser- 
firing delay. The data were background subtracted 
to remove the signal bias resulting from the daylight 
and the detector dark current. The signal profiles 
used in the calculations were limited to the region 
between the range past the full detector turn-on and 
the ground return or the cloud saturation point. 
Because of the difference in pulse lengths of the 
on-line and the off-line laser signals, the smoothing 
intervals of the two signals were adjusted to ensure 
an equivalent 300-m running sampling interval (ver- 
tical range) for DIAL calculations. The effective 
H 2 0 absorption cross-section altitude profile was 
calculated with the Voigt profile representation of the 
line shape, 37 spectroscopic parameters of the H 2 0 
line, and the standard atmospheric model. This 
effective absorption cross-section profile is then used 
to calculate the H 2 0 number density profile by use of 
the DIAL equation |Eq. (2)J. The lidar returns used 
for the DIAL calculations were collected by use of 
shot accumulation times of 1 s onboard the aircraft. 
The 1991 data were further smoothed by use of a 
horizontal average of 6 km ( - 1 min), and for the 


1990 data a variable horizontal interval that ranged 
from 6 km near the aircraft altitude to 30 km near the 
ground was required for the data to retain a good S/N 
over the entire altitude range. The DIAL-derived 
H 2 0 number density profile was converted to a 
mixing ratio profile by use of a midlatitude atmo- 
spheric model, and the error bars in the profile plots 
correspond to the standard error of the average H 2 0 
profile. 

Experimental Data 

In Fig. 8 we show an example of a DIAL H 2 0 mixing 
ratio profile measured off the Virigina coast during a 
flight over the Atlantic Ocean on the night of 15 
March 1990. The stronger H 2 0 absorption line was 
used during this flight, and the range cell applied to 
the return signals was 210 m. This profile repre- 
sents an average of 600 shots, which corresponds to a 
horizontal resolution of approximately 10 km at the 
speed of the Electra. In Fig. 8 the DIAL plot is 
compared with the profiles generated by the two 
dew-point hygrometers as the plane spiraled down 
toward the water. The curves shown in Fig. 8 
indicate very good agreement between the lidar- 
derived H 2 0 distribution and the in situ measure- 
ments, both in the absolute amount and in the 
vertical structure of H 2 0. These data are consistent 
with the high H 2 0 levels and the shallow depth 
(typically < 800 m) expected for a marine mixed layer. 
The higher peak value indicated by the EG&G 300 
hygrometer in the layer located at 600 m may be 
partially the result of an overshoot problem associ- 
ated with the instrument. 

A flight was conducted on 22 March 1990 to 
demonstrate the daytime measurement capabilities 
of the system. As a result of the low levels of 
humidity measured at the WFF ground station that 
morning, the stronger line was again chosen for this 



Fig. 8. Comparison of the DIAL and the DPT (dew-point hygrom- 
eter) H 2 0 volume mixing ratio (parts in 10 3 ) profiles at a location 
appproximately 125 km east of the Virginia coast during a night- 
time flight on 15 March 1990. 
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flight test. Figure 9 presents an afternoon set of 
measurements over Emporia, Virginia. Because of 
the significant variations in the boundary layer, only 
the EG&G 300 hygrometer plot is shown for clarity. 
The difference between the DIAL and the dew-point- 
hygrometer profiles is < 10% over the entire range of 
the measurements. This agreement with measure- 
ments made with in situ probes is typical for the 
DIAL system. These results clearly show the ability 
of the DIAL system to make accurate, high-resolution 
measurements of large variations in the vertical 
distribution of H 2 0. 

After the Electra returned to flight altitude (4 km) 
above Emporia, we flew from Emporia to a point 
approximately 120 km off the coast of Virginia and 
observed the transition from a continental mixed 
layer to a marine mixed layer. Cross sections of the 
aerosol and the H 2 0 distributions encountered on 
this flight path are shown in Figs. 10(a) and 10(b), 
respectively. A 300-m range cell was used, and a 
horizontal average that varied with altitude was 
applied to the H 2 0 data to ensure that the ratio of the 
standard deviation to the mean of the horizontal 
values was < 10%. The aircraft crossed the coast at 
14:43 h local standard time (LST), and the measured 
aerosol and H 2 0 distributions clearly show the de- 
crease in the mixed layer depth from approximately 
1.7 km over land to less than 1 km over water. 
Increased H 2 0 concentrations can be seen near the 
coast in the continental mixed layer as a result of the 
daytime seabreeze circulation, and H 2 0 is seen to be 
positively correlated with aerosol loading. These 
airborne DIAL measurements show the detailed H 2 0 
and aerosol structures that occur in the free tropo- 
sphere and in the mixed layer over different land and 
marine regimes. 

At the end of the 22 March flight test, an overflight 
and spiral descent were performed at NASA WFF 
before landing the aircraft. This procedure provided 



Fig. 9. Comparison of the DIAL and the DPT 300 hygrometer 
mixing ratio profiles over Emporia, Virginia, from an afternoon 
flight on 22 March 1990. 
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Fig. 10. Cross sections of (a) relative aerosol backscatter and (b) 
H 2 0 mass mixing ratio (grams per kilogram) obtained during a 
flight across the coast of Virginia on 22 March 1990. ASL, above 
sea level; LAT, latitude; LON, longitude. 


near-simultaneous H 2 0 profiles from the DIAL sys- 
tem, the EG&G 300 dew-point hygrometer, and a 
radiosonde launched from the WFF ground station. 
Results from these measurements, which are pre- 
sented in Fig. 11, verified the capability of the DIAL 
system to provide an accurate determination of the 
H 2 0 distribution from an aircraft altitude of almost 6 



Fig. 11. Intercomparison of H 2 0 measurements with the DIAL 
system, the DPT 300 hygrometer, and a Vaisala radiosonde over 
the Wallops Flight Facility (WFF) at the end of the flight on 22 
March 1990. 
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km. In addition, the comparison of the radiosonde 
and the hygrometer demonstrates good agreement 
between these two measurement techniques and 
allows us to be confident in using these instruments 
to verify the performance of the DIAL system. 

A daytime flight was conducted across a cold front 
on 2 April 1990 to observe the difference in the H 2 0 
structure on the wet side and on dry side of the front. 
In this case the weaker H 2 0 absorption line was used 
because of the high levels of H 2 0 on the warm side of 
the front. Measurements in the vicinity of the front 
were very difficult to achieve because of clouds and 
storm activity. However, farther from the front, 
clear conditions for DIAL measurements were avail- 
able. One advantage of an airborne system is the 
flexibility of the aircraft in placing the instrument in 
locations that are favorable for DIAL measurements. 
Profiles obtained over Nashville, Tennessee, by use of 
the DIAL system, again with a 300-m range cell and 
the EG&G 300 hygrometer at 16:27 h LST, are 
shown in Fig. 12. The data shown in Fig. 12 are 
representative of the H 2 0 distribution on the dry side 
of the front because the front had passed over Nash- 
ville at approximately 14:00 h LST. It is clearly 
evident from this comparison that the DIAL system 
can make accurate measurements down to very low 
H 2 0-concentration levels, even when using weaker 
absorption lines. It is expected that the accuracy of 
the low concentration level measurements would 
have been greater if the stronger H 2 0 line had been 
used. The H 2 0 cross-section data from this flight 
experiment indicated that the mixed layer height was 
limited to altitudes of 2 km and below on the dry side 
of the front. However, the H 2 0 mixing ratios ap- 
proximately doubled on the warm side of the front, 
and the boundary-layer ceiling increased to approxi- 
mately 3 km. These data show how the DIAL sys- 
tem can measure widely varying H 2 0 distributions 



Fig. 12. Comparison of the DIAL and the DPT 300 hygrometer 
mixing ratios over Nashville, Tennessee, during a daytime flight on 
2 April 1990. The profiles were obtained on the dry side of a cold 
front. 


that extend across different meteorological environ- 
ments. 

In addition to DIAL H 2 0 profiles, the system also 
measures the atmospheric scattering ratio by use of 
the off-line laser return signal, as described above. 
However, because the DIAL system usually operated 
in the nadir mode from an altitude of 4-5 km, it was 
very difficult to find an aerosol-free region for normal- 
ization. As a result, the actual scattering ratios are 
somewhat higher than those that were calculated, 
depending on how clean the atmosphere was at the 
normalization altitude. On the other hand, even 
relative measurements of the scattering ratio can 
provide very useful information. An example of a 
scattering-ratio profile is presented in Fig. 13, along 
with the corresponding H 2 0 plot from a flight on 16 
April 1991. The scattering- ratio scale is delineated 
across the top of Fig. 13. The layer at 2 km, which 
could be a cloud, has a peak scattering ratio of 2. 
A scattering ratio of 2 signifies that 50% of the return 
signal is due to molecular (Rayleigh) scattering and 
50% is backscattered by aerosols ( Mie). Because this 
value could be underestimated, the signal from the 
aerosols is probably greater than that from the 
molecules. The plots shown in Fig. 13 indicate the 
general correlation of elevated H 2 0 levels with re- 
gions of higher aerosol loading. 

The final set of data is shown to illustrate further 
the superior performance and operational capabilities 
of the H 2 0 DIAL system when an APD is used as the 
detector. Figure 14 presents two H 2 0 mixing ratio 
cross sections measured during a nighttime flight on 
18 April 1991. The vertical range cell (210 m) and 
the horizontal averaging interval (2 min ) are the same 
for both cross sections. Figure 14(a) shows measure- 
ments from an altitude of 4.5 km, with the PMT 
mounted in the receiver system. In this case the 
aircraft was flying north along the Eastern Shore of 
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Fig. 13. Comparison of the H 2 0 mixing ratio and scattering ratio 
profiles from measurements performed during the daytime on 16 
April 1991. 


20 September 1994 / Vol. 33, No. 27 / APPLIED OPTICS 6435 




WATER VAPOR (g/kg) 

0 2 4 6 8 10 




(o) MEAN (g/kg) (t>) STANDARD DEVIATION |c) SIGNAL TO NOISE RATIO 

Fig. 15. Performance comparison of APD (aircraft altitude 7.1 
km) and PMT (aircraft altitude 4.5 km) from a flight on 18 April 
1991. The plots show (a) the measured mixing ratio profiles, lb) 
the mixing ratio standard deviation profiles I APD profile adjusted 
for range difference), and (c) the mixing ratio S/N profiles. 
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Fig. 14. Cross sections of H 2 0 mixing ratio distributions obtained 
during a flight along the Eastern Shore of Virginia on 18 April 1991 
by use of (a) a PMT detector in the receiver system (flying north at 
an altitude of 4.5 km) and (b) an APD detector in the receiver 
system (flying south at an altitude of 7.1 km; the APD plot has been 
reversed for easier comparison with the PMT plot). 


Virginia toward a storm system located over WFF 
(the storm edge is at approximately 37.7 °N latitude). 
Of particular interest is the significant increase in 
H 2 0 levels in the vicinity of the storm clouds, which 
has been observed by use of airborne in situ H 2 0 
instruments. 38 For the plot shown in Fig. 14(b), the 
PMT had been replaced with an APD and the plane 
was flying south (the plot has been reversed for easier 
comparison with the PMT cross section) away from 
the storm system at an altitude of 7.1 km. The edge 
of the storm was then located at approximately 36.8 
°N latitude, and the Electra was able to fly above the 
clouds near the edge of the storm. Again, as shown 
in Fig. 14(b), the increase in H 2 0 levels is readily 
apparent in the regions that surround the clouds. 
A 2-min interval ( ~ 12 km) of data from each of the 
two cross-section plots was used to generate a compari- 
son of the performance of the system with the APD 
and the PMT. The approach used to compare the 
H 2 0 distribution measurements is similar to the 
approach used in the receiver detector section, and 
the results are shown in Fig. 15. Although the 
measurements are separated in time by approxi- 
mately 1.5 h, the absolute amount of H 2 0 is relatively 
unchanged. Thus an accurate comparison of the 
performance of the system is possible by use of these 


data. Because the aircraft altitude is higher for the 
APD measurements, the APD standard deviation 
profile was adjusted to compensate for the difference 
in the range from the aircraft to the measurement 
point. It is readily evident from Fig. 15 that the S/N 
for the APD is higher by 50-80% above 1000 m and is 
comparable with the PMT S/N below this altitude. 

Summary and Conclusions 

The NASA LaRC airborne H 2 0 DIAL system was 
designed and developed to perform investigations of a 
wide variety of atmospheric processes by providing 
accurate, high-resolution measurements of H 2 0 and 
aerosol distributions. This paper has described the 
H 2 0 DIAL system and its operation, presented data 
from the initial flight-test experiments, and evaluated 
the system performance. Included in these flights 
were to our knowledge the first airborne mesoscale 
H 2 0 DIAL measurements across a cold front, as well 
as the first measurements of this type in the vicinity 
of a large storm system. Numerous in situ H 2 0 
profile measurements were performed during the 
flight tests for comparison with the H 2 0 profiles that 
were acquired by use of the DIAL system. These 
intercomparisons demonstrate the capability of the 
DIAL system to measure H 2 0 mixing ratios to better 
than ± 10% in diverse atmospheric conditions. Cur- 
rently the system is in a major development phase in 
which the Nd:YAG-pumped dye laser is being re- 
placed with a second alexandrite laser as the off-line 
transmitter. In addition, both alexandrite lasers 
will be injection seeded with diode lasers, and the 
on-line transmitter will be stabilized on an adjustable 
H 2 0 line by use of an active wavelength-modulation 
technique. In the future the H 2 0 DIAL system will 
be used in major field experiments aimed at investiga- 
tion of large-scale tropospheric and lower strato- 
spheric processes, and this instrument will provide 
unique insight into the atmospheric variability of 
H 2 0 and aerosols. 
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and the flight tests of the system: Byron Meadows 
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on interpreting their results; Marta Fenn, Susan 
Kooi, and Greg Nowicki for the DIAL data reduction 
and analysis; Roger Navarro, Dave Pierce, Doug 
Young, and the aircraft crew for their cooperation in 
conducting the flights; Riley Bull and his staff for 
providing the radiosonde data at WFF; and Carl 
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